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ABSTRACT 

Most symbiotic stars consist of a white dwarf accreting material from the wind of a red giant. An 
increasing number of these objects have been found to produce jets. Analysis of archival Chandra 
data of the symbiotic system CH Cygni reveals faint extended emission to the south, aligned with the 
optical and radio jets seen in earlier HST and VLA observations. CH Cygni thus contains only the 
second known white dwarf with an X-ray jet, after R Aquarii. The X-rays from symbiotic-star jets 
appear to be produced when jet material is shock-heated following collision with surrounding gas, as 
with the outflows from some protostellar objects and bipolar planetary nebulae. 
Subject headings: binaries: symbiotic — stars: individual (CH Cygni) — stars: winds, outflows — 
white dwarfs — X-rays: general 



1. INTRODUCTION 

Collimated jets are produced by many accreting as- 
trophysical systems, ranging from proto-stars, through 
stellar-mass black holes in X-ray binaries, to active galac- 
tic nuclei. Symbiotic stars, in which a white dwarf (WD) 
accretes from the wind of a red giant, are probably 
the most recently recognized class of jet-producing ob- 
jects. There are currently at leas t 10-12 symbiotics with 
evidence for bipol ar outflows (Brockso pp et al.l 2004; 
iCorradi et al]l200llh Very little is known about the X- 
ray properties of symbiotic-star jets, as only one symbi- 
otic jet has been re ported in the X-rays (R Aquarii; see 
IKellogg et afll2"f)0l . 

CH Cygni (I = 81?86, b = +15?58) is one of the best- 
studied symbi otic systems ( for a review of CH Cyg's 
properties see iKenvonl 120011) . The multiple photomet- 
ric and radial velocity periodicities present in the system 
sugge st the possible presence of a third star ijHinkle et al.l 
1993). For the distance, we adopt 245 ± 50 pc, from 
the mo de of the Hipparcos p arallax probability distri- 
bution <)Perrvman et al.lll997[) . Evidence for a radio jet 
was first detected between April 1984 and May 1985, 
coincident with a strong radio outb urst and a substan - 
tial decline in the visual magnitude ijTavlor et al.l ll986). 
Subsequent radio and optical observations confirmed the 
exten ded emission as arising from a jet l)Crocker et, alJ 
1200 1|) . and also revealed that the jet position angle on 
the sky changes with time, sugges ting that the jet axi s 
preccsscs with a period of w 6520 d l)Crocker et al.l2 0021. 
Changes in the fast optical flickering in 1997, after a jet 
was produced, have also suggested a direct connection 
betw een the accretion disk and jet production in CH 
Cyg ijSokoloski fc Kenvonl l2003ajh X-ray emission was 
first conclusively detected from CH Cyg by EXOSAT 
on 1985 May 24, at a flux of 1.3 x 10 -11 ergs cm -2 s -1 
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(0.02-2.5 keV: ILeahv fc Tavlorl1 987^ . Subsequent ASCA 
observations revealed a complex X-ray spectrum which 
could be fit eithe r with multiple thermal components 
l)Ezuka et alJl9 98) or a single-temperature the rmal spec- 
trum attenuated by an ionized absorber IWheatlevI 
2001). 

Here we describe analysis of a Chandra observation 
in which we detect, for the first time, the X-ray jet in 
CHCyg. 

2. OBSERVATIONS 

We analyzed an observation of CH Cyg acquired be- 
tween 2001 March 27 6:22 and 20:02 UT, by the Chan- 
dra X-ray Observatory (jWeisskopf et alJ l2002L The ob- 
servation was made in timed/FAINT mode with the high- 
energy transmission grating (HETG) in place, although 
the source flux was too low to obtain a high signal-to- 
noise grating spectrum (the first-order count rate was 
6.3 x 10 -3 counts -1 ). Thus, we did not consider the dis- 
persed spectra for our analysis. The CCDs are sensitive 
to photons in the 0.2-10 keV energy range, and have a 
pixel size of 0'.'492. The zeroth-order image fell on the 
S3 (backside-illuminated) chip of the ACIS-S array. The 
standard data filtering resulted in a net exposure time 
of 47.082 ks. To analyze these archival data, we used 
CIAO 1 software version 3.0.2 and the Chandra calibra- 
tion database version 2.26. 

3. RESULTS 

CH Cyg was detected with a total zeroth order 
count rate of 1.95 x 10 -2 counts -1 . The image 
was not affected by pile- up. The centroid position 
was R.A. = 19 h 24 m 33!076, decl. = +50°14'29'/21 
(J2000.0), which agrees to within O'.'ll with the Hip- 
parcos position of the optical counterpart HD 182917 
ijPerrvman et alJll997f) . This offset is well within Chan- 
dra's 90% pointing accuracy 2 of 0'.'6. The binned 
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Fig. 1. — Chandra & HST images of the jet in CH Cyg. The con- 
tours show the zeroth order Chandra image from the 2001 March 
27 observation, after smoothing with a Gaussian filter of radius 1" . 
The contour values are 0.73, 1.2, 2, 10, 30, 50 and 90 counts per 
pixel. The greyscale shows the HST WFPC2 image from 1999 Au- 
gust 12 in the O III 495.9 nm & 500.7 nm filter (i.e. filter F502N); 
these data are also shown in Fig. 4c of lEvres et al H2002T) . 



zeroth order image of CH Cyg is shown in Fig. ^ 

The source exhibited a faint extension ~ 2-4" to the 
south-southeast. At the source position, 38'.' 1 from the 
optical axis, no significant deviations from the optimal 
point-spread function (PSF) are expected (the 80% en- 
circled energy radius for an on-axis point source is 0'.'685 
at « 1 keV). At the satellite roll angle used, the read- 
out axis was perpendicular to the observed extension. 
To verify the source extension, we generated a sim- 
ulated point-spread function (PSF) with the Chandra 
Ray- Tracer (ChaRT 3 ) in combination with the Model of 
AXAF Response to X-rays (MARX 4 ) simulation soft- 
ware, using the zeroth-order spectrum (see below) as in- 
put. We measured the radial surface brightness profile 
of the simulated zeroth-order image, and the resulting 
point-spread function is shown as a grey ribbon in Fig. [21 
We also measured the observed radial profile of CH Cyg, 
both in the southern region where the source is extended, 
and over the rest of the image excluding this region (see 
Fig. [Q. The extension we observed in the binned zeroth- 
order image (Fig. is apparent as an excess above the 
radial surface brightness profiles of both the simulated 
image and the non-southern directions of the observed 
images. 

We extracted 43 photons from the extended southern 
region, excluding photons within 1'.'72 of the centroid 
position. The probability of observing this many pho- 
tons, given the extent of the simulated point source and 
assuming Poisson statistics, is 1.67 x 10 -13 , correspond- 
ing to a significance of 7.4er. The position angle of the 
extended emission (measured anticlockwise through east 



from north) of the principal component of the photon dis- 
tribution was 164°, consistent with contemporary mea- 
surements of the CH Cyg jet fro m VL A /MERLIN ra - 
dio and HST optical measurements ijCrocker et all 20021 . 
We therefore conclude that the extension observed in the 
Chandra image is due to X-ray emission from the jet in 
CH Cyg. The northern optical jet is too small to be re- 
solved with Chandra, and we see no evidence for extended 
emission in that direction. 

The hardness ratio (calculated from the total counts 
> 2 keV divided by the counts < 2 keV), indicates 
that the extended emission (hardness ratio 0.65) was 
significantly softer than the central source (1.69). The 
photon distribution was, however, qualitatively similar 
to that of the zeroth-order image, with peaks below 
1 keV and around 5 keV (Fig. 0J. To characterize 
the jet spectrum, we fit the data with an absorbed, 
0.18±£j$ keV Raymond-Smith (RS) component for the 
emission < 2 keV and a broad (a « 1 keV) Gaussian at 
5.6^0 5 keV (la errors) for the emission > 2 keV, each 
with a neutral interstellar absorbing c olumn density fixe d 
at 7x 10 20 cm" 2 (from the Hi survey of lStark et alll992|) . 
We found an acceptable fit with a C -statistic (a ppropri- 
ate for low-count per bin regimes; iCashl 119791) of 4.0, 
below which fell only 6% of 1000 Monte-Carlo simulated 
spectra. Extrapolating this model beyond the Chandra 
energy limits, we estimate an unabsorbed bolometric flux 
of 1.2 x 10~ 13 ergscm _2 s _1 , giving a jet luminosity of 
9 x 10 29 ergss" 1 . 

We also fit the total zeroth-order spectrum, taking pho- 
tons from within a circle of radius ss 5" around the 
centroid position. We used photons from an annulus 
5-10" from the centroid position as background. The 
spectrum exhibited two broad peaks, around 1 and 5 
keV, as well as a narrow 6.64 keV line (Fig. Et. simi - 
lar to the earlier ASCA observation (Ezuka et al. 1998). 
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Fig. 2. — Radial surface brightness profiles of the zeroth-order 
Chandra image of CH Cyg. The grey ribbon shows the ±lcr lim- 
its on the simulated PSF from ChaRT/MARX, while the trian- 
gles show the radial profiles calculated using annulae excluding the 
southern extension. The squares show the radial profile for the 
southern extension. Error bars indicate the 1<t uncertainties. 
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Fig. 3. — Chandra spectra from CH Cyg. The zeroth-order spec- 
trum is shown as crosses, with error bars indicating the lcr uncer- 
tainties. The dip in emission between 1—4 keV is not solely a conse- 
quence of a decreased effective area in that region. The histogram 
shows the best-fit model, consisting of 3 Raymond-Smith plasma 
components and a Gaussian to fit the line emission at 6.64 keV. 
The grey histgram shows the jet spectrum, scaled by a factor of 10 
relative to the zeroth-order spectrum (right-hand y-zxis). 

Following those authors, we fit the net spectrum with 
a three-temperature RS plasma model plus a Gaussian 
to describe the 6.64 keV line. We obtained a C-statistic 
of 35.2. 53% of 1000 Monte-Carlo spectral simulations 
gave a lower value, indicating a statistically acceptable 
fit. The spectral parameters were roughly consistent with 
the ASCA values. By extrapolating this plasma-model 
spectrum beyond the Chandra energy limits, we derive 
an intrinsic (neglecting line-of-sight absorption) bolomet- 
ric flux of 1.3 x 10 -11 ergscm _2 s~ 1 , giving a luminosity 
of 9.3 x 10 31 ergss -1 (for a distance of 245 pc). This 
luminosity is more than an order of magnitude smaller 
than typical values from earlier X-ray measurements of 
CH Cyg fe.g. lEzuka et~aDlr99l . 

4. DISCUSSION 

4.1. The X-ray Jet in CH Cygni 

The extended X-ray emission in CH Cyg can be natu- 
rally explained by the collision of a collimated outflow of 
material with surrounding nebular gas. Outflow veloci- 
ties for CH Cyg have been estima ted as 600-2000 km s" 1 
and 1100 kms" 1 in 1 984/1985 llTavlor et all 119861 and 
references therein, and lCrocker et a L||2001| r especti vely), 
and > 1100 kms" 1 in 1998-2000 l)Evres et all 2002). For 
a typical outflow velocity of ~ 1300 kms -1 , the Rankine- 
Hugoniot jump conditions give v s — (4/3)w now ~ 
1700 kms -1 (where v s is the shock velocity and w now 
is the flow speed in the jet) for a strong shock. The 
temperature behind the shock, T ps , would then be 
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where /j, is the mean molecular weight, hp is Boltzmann's 
constant, mp is the proton mass, and we have assumed 
that the gas is ideal. A plasma with a post-shock tem- 
perature of 4 x 10 7 K produces X-rays in the 4-5 keV 
range; although we have insufficient jet X-ray photons to 
perform formal spectral fitting, we detect X-ray-jet pho- 
tons at these energies. The data are thus consistent with 
material heated in a strong shock. 

Softer X-ray photons are also detected from the jet. 
The time scale for hot plasma to adiabatically cool by 
flowing out from behind the shock is given roughly by the 
width of the jet divided by the sound speed. For the jet 
parameters of CH Cyg, this time is a few months. Thus, 
the lower-temperature jet emission (photons around 1 
keV) could be due to material that has expanded per- 
pendicular to the flow of the jet. Alternatively, a range 
of temperatures could arise from multiple shocks in the 
jet. 

Adopting v s — 1700km s _1 and a jet size of 4" 
on 2001 March 27, we can estimate the epoch of 
jet production. For d = 245 pc, we derive a mean 
proper motion for the jet of 1'.'46 yr _1 . The extended 
Chandra emission is therefore consistent with hav- 
ing been produced by material that was ejected in 
mid-1998. At that time, CH Cyg under went a transi- 
tion i nto an optical high state (see e.g. [C^cker^^U 
120011 lEvres et alJ 12002: ISokoloski fc Kenvonl l2003bTT 
The dete ction of ext en ded radio emission in 
1995 by iCrocker et all l|200lD suggests that a 
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Fig. 4. — Comparison of the approximate spectral energy 
distributions for the jet emission from CH Cyg (filled circles) 
and R Aqr (open triangle s). The radio data f or CH Cyg are 
from ICrocker et all I2(XH]) and IKarovsk^e^alJ IHfflSI), a nd for 
R A qr are from iD oughert v et al.l 119951) . ILehto Johnsonj 119921) 
and IKafatos et al l 119891) The UV data for CH Cyg are from 
IK yres et alJ 12002) . and for R Aqr are from lContini fc FormigginJ 
( 2003 ) . The X-ray meas urements for CH Cy g are from this work, 
and for R Aqr are from iKellogg et al.l 120011) . The UV fluxes for 
CH Cyg are not extinction co rrected, and it is not known w hether 
the UV flux for R Aqr from Contini & Formiggini (2003) is ex- 
tinction corrected. The CH Cyg jet X-ray fluxes are absorption 
corrected using the Galactic column only. 
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jet ejection also occurred in association with a transi- 
tion to the optical high s tate in 1992. The other two 
jet ejections, in 1984 /85 l)Tavlor et alJll986|) and 1997 
ijKarovska et al.lll998j) . were associated with transitions 
from high to low optical states. Thus there appear to 
be at least two different scenarios in which jets can 
be produced in CH Cyg. This situation is reminiscent 
of X-ray binaries, which can produce discrete ejections 
in association with X-ray state changes, possibly in re- 
spons e to an extreme physic al change in the accretion 
flow (|Fender fc Hendrvll2000|) . 

4.2. Non-relativistic X-ray Jets 

The only other symbiotic star currently known to have 
an X-ray jet is R Aqr fsee iKellogg et aO 12001^) . Com- 
pared to CH Cyg the R Aqr jet is both more extended 
on the sky (30" for the R Aqr southern jet compared to 
4" for CH Cyg) and physically larger (roughly 6000 AU 
compared to roughly 1000 AU for CH Cyg). Although 
the measured jet fluxes f or the two sources are compara- 
ble. IKellogg et all (f"2001) found that few of the jet X-ray 
photons from R Aqr had energies greater than 1 keV. 
The R Aqr jet spectrum revealed emission lines at 400 
and 550 eV and was generally consistent with a shocked 
plasma out of ionization equilibrium (the plasma temper- 
atures for the northeastern and southwestern jets were 
found to be 2 x 10 7 K and 2 x 10 6 K respectively). A 
follow-up Chandra observation of R Aqr revealed that 
the northern R Aqr X-ray jet is expanding at a veloc- 
ity of sa 600 kms" 1 (E. Kellogg, private communica- 
tion). iContini fc Formigginil l|2003D compiled and mod- 
eled the full spectral energy distribution (SED) of R Aqr, 
and found that the jet emission is due to a combination 
of photoionization and shock heating. In units of vF u , 
where v is radiation frequency and F v is flux density, 
the R Aqr jet emission is strongest in the UV, weak- 
est in the radio, and in between for the X-rays. In Fig. 
[H we plot both the SEP da ta compiled for R Aqr by 
Contini & Formiggini (2003j) and that for CH Cyg ob- 
tained from the literature and our work in the X-rays. 
The jet SEDs for R Aqr and CH Cyg are quite similar, 
despite R Aqr appearing to have a slower shock speed, 
slightly lower post-shock temperatures, and a larger size. 



This similarity suggests that other symbiotic-star jets 
might also be expected to emit X-rays. 

In addition to symbiotic-star jets, two other types of 
non-relativistic jet have been detected in the X-rays: pro- 
tostellar jets (in Herbig-Haro objects) and collimated 
outflows from the central stars or binaries in planetary 
nebulae (PNe). In both cases, the X-rays also appear 
to be due to shock heating as the outflowing material 
collides with surrounding gas. In th e two known X-ray- 
bright protostellar outflows f HH2. | P ravdo et al.| l200ll 
and YL1551 IRS5 = HH 154, iFavata et alJl2002|) . the 
jet X-ray luminosities are a few times 10 29 ergs (similar 
to CH Cyg and R Aqr) and the jet p lasma temperatures 
are ~ 10 6 K. lSoker fc Kastnerl l)2003j) summarized six ob- 
servations of PNe by Chandra and XMM-Newton, four 
of which revealed diffuse X-ray emission with tempera- 
tures of a few time 10 6 K. These four P Ne all show ev- 
idence for high-velocity, collimated flows l|Kastner et alJ 
120031) . ISoker fc Kastnerl (|2003l) suggest that a hidden bi- 
nary companion within each of these PNe accrete from 
the AGB-star wind and produce an X-ray jet. They note 
that the X-ray production mechanism would thus be sim- 
ilar to that in the protostellar outflows mentioned above, 
and the symbiotic star R Aqr. In fact, an additional 
PNe with X-ray-jet emission, Menz el 3, is known to con- 
tain a symbiotic pair at its center llKastner et al l 120031: 
iZhang fc Liill200alSchmeia fc Kimeswengerll20oil . The 
links between the X-ray jets in symbiotic stars, proto- 
stellar outflows, and PNe are thus quite compelling. 
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